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ABSTRACT 

A method is developed to evaluate the magnifications of the images of galaxies 
with lensing potentials stratified on similar concentric ellipses. In a quadruplet 
system, there are two even parity images, two odd parity images, together with 
a de-magnified and usually missing central image. A simple contour integral 
is provided which enables the sums of the magnifications of the even parity or 
the odd parity images or the central image to be separately calculated without 
explicit solution of the lens equation. We find that the sums for pairs of images 
generally vary considerably with the position of the source, while the signed sums 
of the two pairs can be remarkably uniform inside the tangential caustic in the 
absence of naked cusps. 

For a family of models in which the lensing potential is a power-law of the 
elliptic radius, if) oc (a^ + + y'^q^'^Y^'^ , the number of visible images is found as 
a function of flattening g, external shear 7 and core radius a. The magnification 
of the central image depends on the size of the core and the slope (3 of the 
gravitational potential. It grows strongly with the source offset ii (3 > 1, but 
weakly if (3 < 1. For typical source and lens redshifts, the missing central image 
leads to strong constraints; the slope (3 must be < 1 and the core radius a must be 
^ 300 pc. The mass distribution in the lensing galaxy population must be nearly 
cusped, and the cusp must be isothermal or stronger. This is in good accord with 
the cuspy cores seen in high resolution photometry of nearby, massive, early-type 
galaxies, which typically have (3 ~ 0.7 (or surface density falling like distance"^'^) 
outside a break radius of a few hundred parsecs. Cuspy cores by themselves can 
provide the explanation of the missing central images. Dark matter at large radii 
may alter the slope of the projected density; however, provided the slope remains 
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isothermal or steeper and the break radius remains small, then the central image 
remains unobservable. The sensitivity of the radio maps must be increased fifty- 
fold to find the central images in abundance. 

Subject headings: gravitational lensing - galaxies: structure - galaxies: kinemat- 
ics and dynamics - galaxies: halos 

1. Introduction 

For smooth, non-singular lenses, it is well-known that the total number of images is 
odd and that the number of even parity images exceeds the number of odd parity images 
by one (e.g., Burke 1981; Schneider, Ehlers & Falco 1992, chap. 5). When the source lies 
within both the tangential and radial caustics, it is lensed into 5 images. When the source 
lies outside one caustic but inside the other, then it is lensed into 3 images. If the surface 
density at the galaxy centre is cusped and that cusp is stronger than isothermal, then the 
situation changes. There is no radial caustic and there are either 4 or 2 images depending on 
whether the source lies inside or outside the remaining tangential caustic (Evans & Wilkinson 
1998). 

In fact, almost all the known sixty or so gravitational lens systems are 2 or 4 image 
configurations (see Pospieszalka et al. 1999 for details of the gravitational lensing database 
which maintains a hst of candidates). There are two systems in which the presence of a weak, 
central image has been claimed. APM08279+5255 is an ultraluminous broad absorption line 
quasar. Originally, two images were detected serendipitously by Irwin et al. (1998) in the 
optical. Subsequently, Ibata et al. (1999) found convincing evidence for a central third image 
using higher resolution infrared imaging. The source of MG1131+0456 is a radio galaxy. One 
extended radio component is lensed into a ring, another is lensed into two images. There 
seems to be a weak central image at the center of the ring, and this may be the third or fifth 
image of parts of the background radio source (Chen & Hewitt 1993); however, this could 
also be emission from the lensing galaxy. The almost total absence of the central image from 
the known lens systems sets strong constraints on the core radius and the steepness of the 
lensing potential (e.g., Narayan, Blandford & Nityananda 1984, Wallington & Narayan 1993, 
Rusin & Ma 2000). 

The aims of this paper are twofold. On the theoretical side, our aim is to find expressions 
for the sums of the magnifications of any subset of the images produced by multiply lensed 
quasars or galaxies, at least within the framework of a class of flexible and popular models. 
We do this by extending the analysis of our earlier work (Hunter & Evans 2001, henceforth 
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Paper I) to lensing potentials stratified on similar concentric ellipse but this time with cores. 
In doublet (or quadruplet) systems, there are one (or two) even parity images, one (or two) 
odd parity images, together with the missing central image. We provide contour integrals 
which find the sums of the magnifications of the even parity images, the odd parity images 
and the central image separately. On the astrophysical side, our aim is to set constraints 
on the cusp profile and core size of lensing galaxies by requiring that the central image be 
too weak to be detectable. For the 2 and 4 image systems, we investigate the permissible 
core radius as a function of the steepness and the flattening of the lensing potential and the 
external shear. 

The paper is arranged as follows. Section 2 develops the contour integral representation 
needed for lensing potentials that are eUiptically stratified. Section 3 speciahses the analysis 
to models in which the lensing potential is a power-law of the elliptic radius combined with 
external shear of arbitrary orientation. The conditions for triple and quintuple imaging are 
found in terms of the depth and shape of the potential, as well as the position of the source. 
The contour integrals for the sums of the magnifications of the two even parity, the two odd 
parity and the missing central image are evaluated. Section 4 uses Monte Carlo simulations 
to set limits on the lensing potential by requiring the central image to be unobservable. 
Finally, our conclusions are given in Section 5. 

2. Contour Integrals for Observables 

This section develops contour integrals for computing the magnifications of images as- 
suming only that the lensing potential is stratified on similar concentric ellipses. 

2.1. The Lens Equation 

In this paper, we always assume that the lensing potential ^ is stratified on similar 
concentric ellipses with constant axis ratio q with < < 1 so that 

V' = V'(t), T^x' + yV'- (1) 

For a thin lens with potential ^, the lens equation is (e.g., Schneider et al., section 5.1) 

^ ^ X + -fix + -f2y -2x^'{t), 77 = y-F72X - 7iy - — '0'(t), (2) 

where (^,77) are Cartesian coordinates of the source. Here 71 and 72 allow for a constant 
external shear in an arbitrary direction. We always assume that the lens is not circularly 
symmetric (that is, either 5 7^ 1 or 71 or 72 7^ 0). 
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Complex numbers often simplify calculations in lensing theory (e.g., Bourassa, Kan- 
towski & Norton 1973, Bourassa & Kantowski 1975, Witt 1990, Rhie 1997). As in Paper I, 
we shall find it helpful to use 

C = e + i?^= IC|e'^ z^x + \y/q, r ^ zz. (3) 

The lens equation (2) then becomes 

C^Poz + Qz-2z'iI^'{t), (4) 

where 

Po = i[l + + 7i(l - <1% Q = i[l - + 7i(l + g') + 2ig72]. (5) 
It and its conjugate can be written in matrix form as 



where P — Pq — 2iI)'{t), from which it follows that 

z\ 1 f P -Q \ f ( 



(6) 



zj P^-\Q\^\-Q PJ\C 



(7) 



When we form the real quantity t — zz from eqn. (7), we find that the solutions of the lens 
equation also satisfy the real equation 

k{t;CC) = liPC - QC){PC-QC) - [P' - = 0. (8) 

We designate K{t; C, C) = ^ the imaging equation. Its real and positive solutions for r 
provide the image positions z for a given source location Once a solution for r is found, 
the image positions are given by equation (7). Any solution of the original lens equation (2) 
gives a solution of the imaging equation. However, the imaging equation may have negative 
real or complex solutions which do not correspond to true images. 



2.2. Sums of Magnifications 



The lens equation defines a map from the complex {z, 2;)-lens plane to the (C, C)~source 
plane. The Jacobian of this mapping is 
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where the second form has been evaluated using eqn (6). The eigenvalues and eigenvectors 
of J determine the distortion of images. If the eigenvalues are both positive, the image is 
direct; if the eigenvalues are both negative, the image is doubly inverted. These are the even 
parity cases. If one eigenvalue is positive, the other negative, then the image is inverted 
and the parity is odd. The reciprocal of the determinant of the Jacobian of the mapping 
from the real (x, y)-lens plane to the (^, ?7)-source plane corresponds physically to the signed 
magnification of an image (Schneider et al. 1992, chapter 5). With our choice (3) of complex 
coordinates, we find 

db = 

where //j is the absolute value of the magnification and pi is the parity of the image located 
at {xi,yi). For certain positions, detJ vanishes and the magnification is infinite. These are 
the critical points and lines. The caustics arc the images of the critical points and curves 
under the lens mapping (2). Evaluating the determinant of (9) gives 

det J = p2 _ |g|2 ^ (2tP - Qz^ - Qz^)P'{r). (11) 

Differentiating the imaging equation (8) for K partially with respect to r, holding ( and ^ 
fixed, and then using the lens equation (6) to express ( and ( in terms of z and z gives 



9X(t;C,C) -iP'-\Q\ 



dr 



2\ 



[p2 _ |g|2 ^ ^2rP - Qz^ - Qz^)P'{r)] . (12) 



Consequently, we get the following compact expression for the signed magnification of an 
image as 



-g2(p2 _ 


-\Q?) 











Our contour integral representation relies on the special structure of equation (13). Images 
correspond to simple zeros of K{t]C,,C,) and l/{dK/dT) is the residue of l/X at a simple 
zero. We continue the right hand side of equation (13) into the complex r-plane and write 
the sum of the signed magnifications of the images as the contour integral 

Here, C is a contour in the complex r-plane which excludes r = and encloses only the 
simple poles corresponding to whichever visible images we wish to analyse. 



Let us note that our analysis here shows how the methods first developed in Paper I 
can do more than we achieved there. We here exploit them in two new ways. We first 
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relax the restriction that the lensing potential be scale-free and allow it to have a core. 
Cored potentials produce either one, three, or five images depending on the strength of the 
potential, the size of the core and the position of the source. Second, we derive formulae for 
sums of signed magnifications for separate image pairs, and for the central image. Hence, 
the present results are more detailed and informative than those of Paper 1, most of which 
were for sums over the four bright images, consisting of two pairs of opposite parity, and 
which therefore partially cancel. 



3. Power-Law Galaixies with Cores 

We now specialize the analysis to the case of power-law galaxies with a core radius a, 
which have the lensing potential 

^{a^ + rf/^ ifO</3<2, 
^ = { (15) 



A 



ln(a^ + r) if /3 = 0. 



These models were first studied by Blandford & Kochanek (1989) in the context of gravita- 
tional lensing (see also Kassiola & Kovner 1993; Witt 1996; Witt & Mao 1997, 2000; Evans 
& Wilkinson 1998 and Paper 1). They arc the projections of three-dimensional power-law 
galaxies familiar in galactic astronomy and dynamics (Evans 1993, 1994). For example, they 
have been used to model the nearby elliptical galaxy M32 (van der Marel et al. 1994), the 
inner parts of the Galactic bulge (Evans & de Zeeuw 1994), as well as the dark halo of our 
own Galaxy in the interpretation of both the Sagittarius stream and the microlensing results 
(Alcock et al. 1997; Ibata & Lewis 1998). The convergence (or surface density in units of 
the critical density) is 

2?2 (^2 + ^)2-^/2 ■ y^"") 

It is easy to see that 

2^\r) = A{a^ + r)-^"', (17) 

where the parameter B = 2/(2 — f3) and has the range 1 < B < oo. The positive parameter 
A, to which the magnitude of the lensing potential is proportional, can be removed by a 
rescaling. Specifically, we scale all lengths by A^/^. The net effect is to set A — 1 in 
eqs (15)-(17). The dependence of our results on A, which is needed in applications, can be 
recovered by multiplying all powers of a, z and C and their conjugates by the same powers 
of 
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3.1. Images and Caustics 

In this section, we establish the conditions for the numbers and types of images, and 
the forms of the caustics. It is convenient to use the variable t = (a^ + t)^^/^ . It reduces to 
the same variable t used in Paper I in the limit of no core (a ^ 0). When there is a core, 
the physically relevant range of t is restricted to < t < a~^/^. In terms of the variable 

7 = 2^'(r)=t, P^Po-t. (18) 

The imaging equation (8) requires the balance 



1 - aH^ 



p'-\Q\r = (t-hY(t-h)\ (19) 



for those real and positive values of t at which the images occur. The values of ti and t2 
depend only on the flattening and the shear, and are defined by 

ti = Po + \Ql t2 = Po-\Q\. (20) 

It follows from the definitions (5) that Pq > and Pq^ — IQP = g^(l — 7i — 7^) > provided 
that 7i + 7I < 1- Hence, both ti and t2 are positive with ti > t2 and t2 < 1. Fig. 1 
shows how image positions can be found graphically by plotting separately the two sides of 
equation (19). The full curve represents the right hand side of that equation. It depends 
only on the parameters of the lens. The left hand side depends also on the position of the 
source. The three dashed curves display it for three different source positions. The lowest 
dashed curve intersects the full curve curve five times and there are five images. When the 
source is sufficiently close to the center of the galaxy, the left hand side remains small until 
it rises to its asymptote at a~^^^ . A requirement for quintuple imaging is that ti must lie 
to the left of this asymptote, and hence that the core radius must satisfy a < t^^^"^ . Images 
disappear in pairs either when the two curves touch in the (^2, ii) interval, which corresponds 
to a crossing of the tangential caustic, or when they touch in the (ti,a~^/'^) interval, which 
corresponds to a crossing of the radial caustic. If the former happens first, then the three 
images that remain form a "core triplet"; if the latter happens first, then the three images 
form a "naked cusp triplet" (e.g., Kassiola & Kovner 1993). The intermediate dashed curve 
in Fig. 1 represents a naked cusp triplet case of a source which lies outside the radial caustic 
but inside the tangential one. The topmost dashed curve is for a source which lies outside 
both caustics, and gives a single image. There is now just a single image corresponding to a 
single crossing on (0, ^2)- The three different source positions are shown in Fig. 2a. 

The maximum number of images diminishes as the core radius increases. If the core 
radius satisfies a < t^^^^, then quintuple imaging is possible, and there are two caustics. 
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For core radii in the range < a < ^2 ■> the vertical asymptote of the dashed curve in 

Fig. 1 hes between ti and t2- There are then three images for the source sufficiently close 
to the center of the galaxy that it lies within the tangential caustic, and there is no radial 
caustic. For still larger core radii satisfying a > then only a single image can occur and 

there are no caustics. Appendix A justifies the statements concerning images and caustics, 
and gives equations for determining the caustics, and an approximate formula for the radial 
caustic for small a and B <2. 

Image positions and points on caustics must generally be determined numerically. Cusps 
arc an exception because they occur in pairs when either t = t2 or t = ti is a triple root of the 

D /O 

imaging equation. The sole caustic is tangential and of lips type for the range < a < 

t2^^^ as in Fig. 2c because it has only the t = t2 pair of cusps. When a < ti^^^ and quintuple 
imaging is possible, there are two pairs of cusps and three possible configurations of caustics. 
Either both pairs of cusps lie on the tangential caustic, which is of astroidal shape, or else 
both caustics are lips-shaped with the radial caustic lying within the tangential and oriented 
oppositely to it (Schneider et al. 1992, Section 8.6.1). This double hps configuration occurs 
when 



as in Fig. 2a. The t — t2 pair of cusps are then naked. The stage at which those cusps become 
naked is described by the coupled pair of equations ( A5) , which must be solved numerically. 
The results are plotted in Fig. 3. The tangential caustic has a pair of naked cusps in the 
regions of parameter space below the curve for the appropriate value of B. Above that curve, 
the tangential caustic has four cusps and lies wholly within the radial caustic as in Fig. 2b. 
Below that curve, but where t2/ti is sufficiently large that the condition (21) is violated, the 
two caustics intersect as in Fig. 2d, and two of the four cusps of the tangential caustic are 
naked. 

We note that our results extend previous calculations in the literature. For example, 
Kassiola & Kovncr (1993) give the conditions for multiple imaging in the isothermal case 
(/3 = 1) in the absence of shear (71 = = 72), while Evans & Wilkinson (1998) give the 
results for the scale-free cases (a = 0) with on-axis shear only (72 = 0). 




(21) 



3.2. Magnifications of Pairs of Images 



Changing the variable in the contour integral (14) for the magnification to t gives 




(22) 
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Fig. 1.— Graphs of i^2(0/(l - «^^^) (dashed curves) and -Ki{t) ^ (t - ti)^{t - ^2)^ {Ml 
curve) for the lens with B — 1.5, q — 0.8, a — 0.84, 71 = 0.1, and 72 = 0. The intersections 
of the curves of different type correspond to images. The three dashed curves, from the 
lowest up, are for the three source positions ( = \(\ exp(47ri/9) for \(\ = .08, \(\ = .15, and 
\C\ = .225. These source positions give five, three, and one images respectively. Only for such 
a large core radius can a figure be drawn in which all the intersections and the asymptote 
t = a^^/^ are visible on the same scale. 
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Fig. 2. — The four different caustic configurations wliicli can occur with the cored elliptical 
potential (f5). Case (a) is for the double lips caustic for the lensing potential of Fig. 1. The 
dots show the three different source positions, aligned and with increasing \(\, for the dashed 
curves of Fig. 1. Case (b) is for the B = 2, q = 0.8, a = 0.1, 71 = 72 = case for which Fig. 4 
is plotted. The single hps case (c) differs from (a) only in having the larger core a = 1.0 for 
which there is no radial caustic. Case (d) differs from (b) only in having the smaller axis 
ratio q — 0.7. This model lies below the B — 2 curve in Fig. 3 and has naked cusps. 
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Fig. 3. — The regions of parameter space in which the tangential caustic has naked cusps 
are those which he below the curves plotted. They grow as B increases. The values of t\ 
and ^2 depend on the shear. The approximations in the labels are those for negligible shear. 
The radial caustic balloons in size as a — > for S < 2, as predicted by equation (A7). 



We write the denominator term as 



(1 - aH'')K = i^i + K2, 
where we have labelled its two parts as 

K, = -{l-aH^){P'-\Q\r = -{^-aH''){t-Unt-hr, 

K2 = t^'iPC-QOiPC-QO- 



(23) 



(24) 



The component K2 is the same as in Paper I, whereas Ki differs from the Ki — — (t — ti)^(t — 
^2)^ of Paper I only through its extra (l — a^t^) factor, which is unity when there is no core. 

In Paper I, wc were interested primarily in sums over the four bright images. However, 
the contour integral (22), like the corresponding result in equation (23) of Paper I, applies 
equally well to any subset of the images provided that the contour C is chosen to enclose 
only the zeros of (1 — aH^)K which correspond to that subset. If a < t^^^'^ for instance, we 
define a contour C2 to be a simple closed contour in the complex t-plane which encircles the 
two direct images of even parity. These approach t = t2 as |C| ^ 0, and so the contour C2 
also encloses t2 which lies in between, but it excludes all the other images, as well as t — 
t = ti, t = a~^/^ and any complex zeros of (1 — o^t^). Then, we obtain the formula 



2 direct 
images 



2 direct 
images 



2m 



dt 



C2t{t-h){t-t2){l-aH^) (l + f)' 



(25) 



We require that C2 keep a finite distance from ti, t2, and every zero of (1 — aH^). Then, 
with \t\ bounded above, and |t — |i — ^2], and |1 — aH^\ bounded from below, it follows 
that 



t 



B 



< 



1 - aH^ 



2+B 



aH 



{t-t,y 

\t-h\' 



{t - hY 
12 " 



+ 



t-U 



(26) 



and hence \K2/K\\ < 1 on C2 for sufficiently small \C,\. We have here used the result 



{PC-QO{PC-QO = [\{t2-t)+v{h-t)][\{t2-t)+v{t^-t)] = \\Wt2-tf+\vWt^-t)\ (27) 



where A and u are the following two linear combinations of the complex source coordinates 
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which have the property |Ap + |i/p = With \K2/Ki\ < 1, we can now expand the 

(1 + K2/Ki)~^ term in the integral (25) as a geometric series to get 

-q'^B r dt 



2 direct 27ri J^^ t{t - h){t - t2){l - aH^) 

images 



q'B r t^^-'[mt-t2r + w\\t-t,ry 



j^i JC2 (t - ti^'+'it - ^2)2^+^1 - aH^y+^ 

The restrictions that we imposed earlier on C2 guarantee that this series converges for suf- 
ficiently small Id- Every integral in it can be evaluated from its residue at t = ^2 alone 
because that is the only singularity within the contour €2- The leading order term is the 
residue of the simple pole at i = ^2, which is positive because both parities are even. The 
full expansion for the sum of the two magnifications is 

q^B 

^ ^-^ 2|g|t2(l-a^tf) +^^^^^'^^'"''^'"^- 

2 direct 1^1 2 / 

images 

Here, C2 is a correction term which vanishes if the source is exactly aligned with the center 
of the lensing galaxy. More generally, it takes the form: 

C2{h,t2, a^; \'^) = q'Bj2J2[rn) IM'^'^l^'^'^'Mti, ^2, a'-j, 2m, 2j - 2m), (31) 

j=l m=0 ^ ^ 

where we have defined 

1 / t^^~^dt 
R2{t,,t2,a^-JXk) = -—f^^ it-t,Y^^it-t2f^\l-aHBy^^ - ^^'^ 

In fact, we can obtain closed form expressions for all of these integrals without any further 
residue calculus by simple partial differentiation because 

1 d'^ 

R2{ti,t2,a^;j,£,k) = -—R^{ti,t2,a^;j,£,0) (33) 



1 d'' 



^2 



If the conditions for quintuple imaging are satisfied, then there is a pair of inverted 
images which approach i = as |C| — > 0. We can perform a similar analysis for a simple 
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closed contour Ci in the complex t-plane which encloses the images and ti, but not t — 0, 
t — t2, t — a~^/^ or any complex zero of (1 — aH^). The result is that 



—qB 

E E ^^ = M^7T3^ + ^i(^^'^^'«';^'^)' (34) 

inverted 2 inverted i i v i / 

images images 



where the correction term Ci is 

oo J 



C,(h,t2, a^; X,u)^q^Bj2J2{} \M'''"H''~^"'Ri(ti, ^2, a^'J, 2m, 2j - 2m). (35) 



j=l m=0 

Here, we have defined 

,2. 



Ri(ti,t2,a ■,j,£,k) = 



1 r t^i-^dt 



1 

-—R^{ti,t2,a^;j,0,k) (36) 

ami 



1 tf-^ 



(ti - t2)^+^ (1 - a^tf )i 
= R2it2,ti,a'^;j,k,i). 

In fact, the coefficients Ri{ti,t2, a^; j, ^, A;) and 7^2(^1,^2, a^; j, £, /c) can be expressed as finite 
sums, as is demonstrated in Appendix B. 

Neither series (30) nor (34) remains convergent at the tangential caustic where one of 
the direct and one of the inverted images merge; the magnifications of the merging images 
become infinite and one cannot then construct the Ci and C2 contours. However, the sum 
of series (30) and (34) can remain convergent at the tangential caustic because the infinities 
of the signed magnifications of the two merging images cancel. The series (34) ceases to 
converge at the radial caustic where an inverted and a doubly inverted image merge, but the 
series (30) for the direct images converges there if, as in Fig. 2a, it lies inside the tangential 
caustic. 

Let us note some interesting special cases of the preceding formulae. In the scale-free 
limit (a = 0) when B = 1 or B = 2, then remarkably the correction terms have the property 
that Ci — C2 = 0. As first discovered by Witt & Mao (2000; see also Paper 1), the sum 
of the four signed magnifications is then an invariant completely independent of the source 
position 

4 images 
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Even though this result is not exact for other values oi B, it is often an extremely good 
approximation. If the core radius is non-zero, then the sum of the signed magnifications 
becomes 



E 

4 images 



B 



a 



+5+1 ^ 
^2 J 



tl-t2 



(38) 



and this approximation remains remarkably accurate inside the radial caustic. An example 
is the cored isothermal lens of Fig. 2b. The sum of the four signed magnifications is constant 
to within tenths of one percent over all of the region inside the tangential caustic except very 
close to the cusps. However, the smallness of the dependence of that sum on the source posi- 
tion is due to a near-cancellation of the correction terms in the two image sums. Fig. 4 shows 
how much the sum of the magnifications of the two direct images varies over the same region. 
Numerically, the leading coefficients in expansion (31) for C2 are i?2(ti, ^2, o^'-, 1, 2, 0) = 13.77 
and i?2(^i! ^2; Q'^; li0,2) = 21.64. The sums -Ri + R2 for the same sets of indices are -0.01. 
Though numerical values vary, Ri + R2 is typically at least two orders of magnitude less than 
R2 for lenses with small cores and with an inner tangential caustic. In Section B.3, we show 
how near-cancellation can occur more generally from the contour integral formula (14), and 
is not a peculiarity of power-law galaxies. However, this behavior is not universal. The sum 
of the four signed magnifications has large variations over the region inside the inner radial 
caustic of Fig. 2a, while the sum of the magnifications of the two direct images changes little 
in the middle third of that region. That middle third lies well inside the tangential caustic 
at which the direct image sum becomes large. 

The separate sums of the pairs contain more information than their signed combination. 
Physically speaking, it is the total magnification (that is, the sum of the magnifications) 
which is more interesting, but this varies considerably with source position. The minimum 
magnification Amin in the scale-free limit (a = 0) is 



4 images 



D 



Pn 



i-ii-ii\Q\ 



(39) 



This corresponds to the case when the source is aligned with the center of the lensing galaxy. 
It is a good approximation only when the source offset is small. For non-zero core radius, it 
becomes 



A. 



4 images 
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1 - aHf 



1 - aH^ 



(40) 



In the next section, we discuss another nearly constant four-image sum which reduces to the 
sum of the unsigned magnifications when the source is aligned with the center of the lensing 
galaxy. 
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25 




Fig. 4. — The surface given by the sum of the magnifications of the two direct images for 
B — 2, q = 0.8, a — 0.1, 71 = 72 = over the first quadrant of the (^,ry)-plane inside the 
tangential caustic. That caustic is drawn as the dashed curve on the base of the plot. The 
full curves are the intersection of the surface with the planes r) — 0, and evenly spaced planes 
^ = const. The filled circles denote the points at which 77 is 0.96 of its value at the caustic. 
Curves which do not end with filled circles exceed 25 before getting that close to the caustic. 
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3.3. Configuration Moments 

The contour integral method can also be used to calculate configuration moments, that 
is, sums over the images of products of the signed magnifications with position coordinates. 
Those sums are obtained by adding the complex position coordinates from (7) to the sum 
(14) to obtain 

We expand for moments in the same manner as for the magnifications, with the result that 



2m j^Jc. [{t-h){t-h)f^^^+-+\l-aH^y+^ 

images 

(42) 

After the two numerator factors are expanded binomially, we see that every integral is of the 
same form as (32), and so, provided m + j > and n + j > 0, we obtain 

oo 2j+m+n 

E ^i^zTz'i^{-lr+^q^BY, E ^2(ii,i2,a';j,^,2j+m + n-£) 

2 direct j=0 i=0 

images 

min(€,m+j) , i A / i A 



X 

k=max[0,£—n—j) 



The corresponding sum for the indirect pair is 

.oo 2j+m+n 

E l^izTz^^i-ir^'^+VBY, E Ri{ti,t2,a';j,i,2j + m + n-i) 

2 inverted j=0 i=0 

images 



min(^,m+j) 

X 

fc=max(0,£— n— j) 



In either case the same R\ and R2 functions as arise with the magnifications, and which are 
evaluated in Appendix B, are all that is needed. 

Whereas the present paper is mostly concerned with magnifications, our modeling in 
Paper I made use of configuration moments too. The analysis of this paper, and especially 
that of Appendix Section B.3, shows that four-image sums of moments may also be far 
more uniform than the separate sums (43) and (44) because of cancellations. However it 
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is possible to mitigate the cancellations which arise from signed magnifications by using 
specially contrived configuration moments. As an example, consider the zeroth order m = 1, 
n = —1 moment [see eqn. (48) of Paper I] in which the magnifications are weighted by the 
complex exponential e^**, where z = |2;|e**. The angle $ for each image is arctan(i/j/ga;j)This 
factor exactly compensates for the different signs of the magnifications when the source offset 
is zero and the two kinds of images lie on perpendicular lines through O. We obtain the 
j — terms directly from the integral (42) and its pair because the binomial expansions 
used to obtain equations (44) and (43) are not vahd when n + j < 0. This gives us two 
R{ti,t2, a^; 0, 0, 0) integrals whose sum gives 

1 - aH^ ^ l-aH^\ ' ^^^^ 

because \/ \ — —vj v — Ql\Q\. The sum (45) is minus the sum in equation (40) when 
there is no off-axis shear and Q is real because e^'* = 1 on the major x-axis where the 
inverted images initially lie and e^** = — 1 on the minor y-axis. Off-axis shear rotates the 
configuration and Ql\Q\ is then the complex factor which makes ^ A*iPi(Qe^**/|(5|) real. 

4 images 

This four-image sum can remain nearly constant as the source offset changes because the 
changing angular positions of the images counteracts their changing magnitudes. For the 
cored isothermal lens for which Fig. 2b and Fig. 4 are plotted for example, the m = 1, 
n = —1 moment of equation (45) displays the same near constancy over the region inside 
the tangential caustic as does the sum of the four signed magnifications. Estimates of q and 
the orientation of the lens must be combined with observed positions to evaluate the sum in 
(45), but the result should be real when multiphed by Q/\Q\. 
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3.4. The Magnification of the Central Image 

There is no fifth central image in the absence of a core for B <2. When there is a core, 
the central image is given by the root of the imaging equation (23) for which t a^"^/^ and 
T ^ as Id — 0. We can obtain formulae for its signed magnification involving another set 
of it!-functions, defined by 

Kit,, i, A.) = — _,2,B).+i (46) 

where Ca is a loop enclosing t — a~'^/^ but not t-i or t2 or or any other zeros of (1 — aH^). 
These functions can be calculated by evaluating residues aXt — a"'^!^ . This is a simple pole 
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when J = 0, and we then get 
The signed magnification of the central image is given by 



(47) 



(48) 



[l-a2/Sti][l-a2/%] 
where the correction term Ca(^i, ^2,0^; A, i/) has the form 

00 j / ■ \ 

CaihM, a'; A, I/) = 5^ ^ r I A|2-|i.|2^-2™i?„(ti, t2, a"- j, 2m, 2j - 2m). (49) 

j=l m=0 ^ ' 

The lowest order term in (48) is positive if a^"^!^ > ti when this image is doubly inverted, 
negative if ti > a'"^/^ > t2 and the image is simply inverted, and positive if ^2 > and 
it is the only and direct image. 

When a is small, as in our applications, so that a"^/^ ^ ti and the central image occurs 
for a large value of t, equations (48) and (49) can be approximated by 



l^aP = l^a = q^a^'"" + q^B ^ |C|'^i?a(0, 0, a^; j, 2j, 0). 



(50) 



These special cases of integrals (46) can be evaluated by residues (change to a^t^ — 1 as 
integration variable) and give a series expansion for the magnification of the central image 
as 

00 

l + j:a,(|C|V(-^)/^)^ 



1 
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r2(j+i) .1 

B 






'2(j+l) , 1 


- j 



(51) 



The series is especially simple for B = 2 when it is binomial and gives 
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There are also two cases for which the series is hypergeometric and can be summed explicitly. 
They are B = 1, for which 

4?2a4 



and B — 4:, for which 



fia = q a 



ICP^ l + |ClV2a^ 



a Vl + lClV4a2 



(53) 



(54) 



The reason for the simphfications is that they are most easily obtained directly from the 
basic magnification formula (22) after its integrand is approximated for large t to give 



q'^B r dt 



2m %jS[tB-2l(^\2-i + aH^y 



(55) 



The root for the central image and its residue can be calculated directly without any need 
for series expansion for the three special values of B. More generally, series expansion is 
needed. This analysis makes it clear that the basic requirement for the approximation (51) 
to be valid is that t be large at the central image. 



for B <2 because all its coefficients aj are then positive. Numerical evaluation, as in Fig. 5, 
shows that Ha remains an increasing function for larger B. Those plots extend to 96% of the 
radius of convergence. The series converge for 



as follows from using Stirling's formula to approximate the gamma functions and standard 
convergence tests. Hence the range of usefulness of the series depends considerably on the 
value of B. For B < 2 and a small, the radius of convergence is large, and the series converges 
out to the radial caustic as approximated by equation (A7). The magnification of course 
becomes large as the caustic, at |C| = 1 for i? = 2 and \(\ = 1 /2a for B = 1, is approached. 
For i? > 2 on the other hand, the radius of convergence is caused by the breakdown of the 
approximation that the central image occurs for large t and, as with equation (54), /i^ remains 
finite as the condition (56) is violated. That radius of convergence in \(\ is proportional to 
the positive power a^~^/^. Consequently, the approximation (51) is useful for only a small 
part of the region within the radial caustic when the core radius a is small, and /^a, which 
is growing with source offset on a short length-scale, can continue to become considerably 
larger before the radial caustic is approached. 

Numerical tests shows that the approximation (51) works well where the theory predicts 
that it will; over much of the region within the radial caustic for B <2, but only in limited 
central regions for i? > 2. It underestimates Ha because of the neglect of denominator terms 
such as those at the leading order in equation (48). That underestimate can be by 10 % or 
more for a = 0.1, but accuracy increases substantially for smaller a. 



It is evident that the series 



(51) shows fia to be an increasing function of ICP^^''^ ^^^^ 




(56) 
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Fig. 5. — Sums of the series (51) for the scaled magnification of the central image as functions 
of its scaled argument out to 96% of their radii of convergence. Those radii vary with B as 
in equation (56). The location of the radius of convergence approximates the radial caustic 
for the top two curves with B < 2. The near-linear growth with increasing source offset of 
the lower two curves for S > 2 is a little less than, but growing towards, the B > 2 estimate 
|^|4/(s-2) equation (51) of Paper I. 
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Fig. 6. — The percentage of quintuplets with a visible central image is plotted against the 
slope of the surface density. If the model is nearly singular, then 2/5 > 1 for cusps steeper 
than isothermal, while 2/B <1 for cusps shallower than isothermal. The four panels show 
results for different values of the core radius a. There is no fifth image when a = and 
2/5 > 1. In each panel, the full line denotes q = 0.7, the dotted line q = 0.8 and the dashed 
line q = 0.9 models. The threshold is 1%. 
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Fig. 7. — The percentage of core triplets with a visible central image is plotted against the 
slope of the surface density. The four panels show results for different values of the core 
radius a. In each panel, the full line denotes q = 0.7, the dotted line q = 0.8 and the dashed 
line q = 0.9 models. The threshold is 1%. 
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Fig. 8. — The percentage of naked cusp triplets with a visible central image is plotted against 
the slope of the surface density. The two panels show results for different values of the core 
radius a. In each panel, the full line denotes q — 0.7, the dotted line q — 0.8 and the dashed 
line q — 0.9 models. The threshold is 1%. 
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4. Application: The Missing Central Images 

Here, we study the problem of missing central images from the standpoint of the struc- 
ture of the lensing galaxy. In optical lenses, the experimental constraint is weak, as any 
central image can be masked by emission from the lensing galaxy. Radio-loud lenses provide 
much stronger constraints, as they have typically been probed with high dynamic range 
radio maps. The Cosmic Lens All-Sky Survey (CLASS) is the largest statistically homoge- 
neous search for gravitational lenses (e.g., Myers 1999). The survey sample contains ~ 10^ 
flat-spectrum radio sources. In the best cases, such as B0218-I-357 (Biggs et al. 1999) and 
B1030-I-074 (Xanthopoulos et al. 1998), the magnification ratio of the faintest to the bright- 
est image r is constrained to < 0.1% from the absence of a detectable central image in the 
map. As a typical detection limit in the below calculations, we adopt r ~ 1% for central 
images. This is appropriate for radio lenses, but not for optical lenses. 

We perform Monte Carlo simulations using the power-law galaxies (15). Sources are 
placed randomly within the outermost caustic for choices of q and core radius a. The 
imaging equation (19) is solved numerically to find the roots t of the images. We then 
evaluate the magnifications using (13), and work out the ratio of the brightest to the faintest 
image. Repeating this many times gives us the raw probabilities of observing a central image. 

Figs 6-8 show the probability of observing a central image as a function of the slope of the 
projected density (2/5 = 2—/?) for the case of the quintuplets, the core triplets and the naked 
cusp triplets respectively. The panels show how the raw probabilities depend on the core 
radius a, while the different lines in each panel show different flattenings. The models have 
a singular density profile if a = 0, and are very nearly singular when a is small. We describe 
the parameter 2/B as the slope index, as it controls the fall-off of the projected density. For 
the isothermal cusp, 2/B = 1 and so the projected density k ~ distance"^. Models steeper 
than isothermal have 2/B > 1, models shallower than isothermal have 2/B < 1. 

There are a number of things to notice in Fig. 6. As the core radius a is diminished, the 
regime in which the central image is visible shrinks. Visible central images occur only for B > 

2 when the central image is stronger as the simple estimate fia ~ q^a^^^ from §3.4 predicts. 
The full curves for the fiattened (g = .7) models lie above the dashed curves for the rounder 
{q = .9) models. This is due to a stronger dimming of the brightest image with increasing 
fiattening which outpaces the dimming of the central image. All magnifications share a g^- 
dependence which comes from the basic magnification formula (22). The brightest image is 
direct for vanishing shear and zero source offset, and has magnification //i ^ q'^B /{A\Q\t2) = 
B /2{1 — q^). Hence we estimate the ratio of the central image to the brightest image at zero 
source offset to be ^ 2q^{l-q^)Ba^/^ . This ratio increases with increasing flattening. 

For small a, the magniflcation of the central image does not vary greatly with increasing 
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source offset witfiin tfie quintuplet region for S < 2. This region inside the tangential 
caustic is but a small part of that within the radial caustic, and hence the abscissa of Fig. 5 
varies little. Conversely, the brightest image strengthens considerably with increasing source 
offset as it gets closer to the tangential caustic. The result is that the ratio /i^/ /Xi of central 
to brightest magnifications decreases with increasing source offset for B < 2. For > 2 on 
the other hand, Ha grows so strongly across the quintuplet region that it can outpace the 
growth of /Xi. The net result is that the flatter the potential, and the larger the core radius, 
then the greater is the likehhood that the central image is bright enough to be visible. 

Fig. 7 shows the raw probability of observing a central image for core triplets. The 
behavior of the curves as a function of slope index 2/B has the following explanation. If 
2/ B is too large, then the central image is highly demagnified and so the raw probability is 
vanishingly small. As the slope index 2/B diminishes, the central image becomes brighter 
and the probability rises quickly to 100%. The brightest image is generally the remaining 
direct one. It dims with increasing source offset whereas the central image brightens. The 
tangential caustic also grows in size as q decreases for constant slope index. That is the 
main reason why central images are more visible with increased fiattening for B > 2; the 
average source offset in the core triplet region between the tangential and radial caustics 
is then larger, and hence is significantly larger. As 2/B decreases further, the astroidal 
tangential caustic grows, and the area in the source plane generating core triplets diminishes 
and eventually vanishes, as the tangential caustic becomes larger than the radial caustic as 
in Fig. 2c. The smaller q is, the sooner this happens and, as Fig. 7 shows, the size of the 
core has little effect on the stage at which the core triplet region disappears. The second 
and third terms of the inequality (21) become equal at that stage and, in the absence of 
shear, 2/B is then (1 — q'^)/{l — a^) (See Appendix A). If the core radius a > 0.01, then 
even models steeper than isothermal {f3 < 1) can provide observable central images. For 
smaller core radii (a < .01), the visible central images are confined to models less steep than 
isothermal. Notice that the constraints on the maximum possible steepness of the lensing 
potential provided by the missing image are stronger for doublets than quadruplets once 
a ^ .01. For small core radii, the fifth image in quintuplet systems is significantly more 
demagnified than the third image in triplet systems. However, for larger core radii, it is the 
quadruplets that provide the stronger restriction. 

Fig. 8 shows the raw probability of observing a central image for naked cusp triplets. 
Naked cusps are much less abundant for elliptic potentials as opposed to elliptic densities 
(Kassiola & Kovner 1993). Almost as soon as naked cusps appear, all three images are 
of roughly similar brightness and they are all detectable. Hence, the raw probability of 
observing a central image shows a swift transition from nearly 0% to nearly 100% as soon 
as naked cups become possible. The Monte Carlo results are consistent with the transitions 
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given in Fig. 3. Equations (A5) predict that the transitions occur at the values 2/B — 
0.565, 0.913, 1.157, for a = 0.1 and 2/B = 0.542, 0.850, 1.004, for a = 10"^ There are no 
strong observational candidates for naked cusp triplets, and so we must conclude that Fig. 8 
sets a firm lower limit on the slope index 2/B. This must be larger than the critical value 
which permits naked cusps (given in Fig. 3), otherwise naked cusps would be common. 

To compare with data from surveys, we must allow for the amplification bias. Lens 
systems with a high total magnification ^ are preferentially included in a fiux-limited survey 
(e.g., Turner 1980, Turner, Ostriker & Gott 1984). The fiux distribution of the sources in 
CLASS is well described by dN/dS oc S*"^'^, where N is the number of sources with fiux 
greater than S (Rusin & Tegmark 2001). For a fiux limited sample, the probabilities that 
take into account amplification bias are 



where A denotes the area enclosed by the caustics in the source plane for which the central 
image passes the threshold (e.g., Rusin & Ma 2001). Fig. 9 is the analogue of Fig. 7 when 
amplification bias is taken into account. Only the rightmost branch of the curve is plotted, 
as this is the most relevant for constraining the core size and the slope index. Notice that 
the effects of the amplification bias cause only slight changes in the shapes of the curves. 
Once the core radius falls below a ^ 0.01, then irrespective of the fiattcning, the lensing 
potential is constrained to be at least as steep as isothermal to ensure that the fraction of 
triplets with an observable central image remains low. As Fig. 10 shows, this conclusion 
remains vahd even in the presence of shear. Shear has httle direct effect on the magnitude of 
the central image as equation (51), which is independent of shear, predicts. However, shears 
of the order of 0.1 cause the inner tangential caustic to grow significantly in size for the 
2/B values at which the curves of Fig. 10 rise sharply, while the radial caustic merely tilts 
a little. The greater visibility of the central image is again because the average source offset 
in the diminished core triplet region has become larger, now as a result of shear. The total 
shear produced by internal misalignments, large-scale structure and neighboring galaxies is 
typically constrained by I7I ^0.3. 

What restriction is implied on the core radius a from the data on radio lenses? The 
CLASS survey found no triples, but 7 doublets. There is perhaps only 1 definite triple 
system (APM08279+5255) out of a total of ~ 50 doublets and triplets on Pospieszalka et 
al.'s gravitational lensing database (which contains both radio and optical lenses). So, the 
probability of detecting a triple with an observable central image is certainly very low. In 
this paper, we take it to be ^ 2% for radio lenses (i.e., at the 1% threshold). Using Fig. 9, the 
top left panel showing curves for a = 0.1 suggests that detectable central images are common 
for isothermal (or nearly so) galaxies. The smallest core radius a seemingly compatible with 
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Fig. 9. — The percentage of triplets with a visible central image is plotted against the slope 
of the surface density. This figiu^c incorporates the amplification bias, i.e., the tendency of 
high magnification configurations to be preferentially included in a fiux-limited sample. The 
threshold is 1%. The increase of visibility with increasing flattening is now due primarily to 
the increasing size of the inner tangential caustic, which subtracts an area of weaker central 
images from the triplet region. 
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Fig. 10. — This figure shows the effects of external shear. For models with a — 0.01, q = 0.8, 
the percentage of triplets with a visible central image is shown for a range of values of the 
shear components. The threshold is 1%. 
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Fig. 11. — The percentage of triplets with a visible central image is plotted against the 
threshold (expressed as a percentage). The effects of amplification bias are included. The 
four panels show results for different slope indices (0.75, 1.0, 1.25 and 1.5). In each panel, 
the full line denotes q — 0.7, the dotted fine q — 0.8 and the dashed line q — 0.9 models. 
The core radius a = 0.01. 
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the missing central images is a 0.01 if the galaxy is isothermal. At a typical lens redshift 

of Z\ ~ 0.5 and a typical source redshift Zs ~ 2.5, this corresponds to a physical size of 
~ 50 pc (using formula (11) of Young et al. (1980) and formula (2.4) of Wallington & 
Narayan (1993)). The result holds good for a flat, matter-dominated Friedman-Robertson- 
Walker universe with a Hubble constant Hq = 75 kms^^Mpc^^. Suppose instead we use the 
currently popular cosmological model in which ~ 70% of the mean energy density required to 
make space flat is in the form of material for which gravity acts repulsively and the remaining 
~ 30% is carried by coUisionless massive particles of some type. Then the physical size of 
the core radii of isothermal galaxies is fa 70 pc. These results are good for 2/B — 1. Prom 
the top right panel of Fig. 9, the dimensionless core radius can be increased to a 0.05 if the 
galaxy has a larger slope index (say 2/B ^ 1.3, which we will argue shortly is appropriate 
for giant ellipticals). In the same two cosmological models and with the same assumptions 
as to typical source and lens redshifts, this gives physical sizes of the core radius of ^ 260 
pc and ~ 350 pc respectively, consistent with almost all the central images being absent. 

Most of the optical depth to strong lensing resides in giant elliptical galaxies. Faber et 
al. (1997) analyse Hubble Space Telescope surface photometry of nearby ellipticals. They 
provide convincing evidence that giant elliptical galaxies have cuspy cores with steep outer 
power-law profiles and shallow inner profiles separated by a break radius. This is in contrast 
to low luminosity ellipticals, which have power-law surface brightness profiles. There are 26 
giant ellipticals with cuspy cores presented in Table 2 of Faber et al. The mean value of the 
asymptotic outer slope of the surface brightness is 1.28 with a standard deviation of 0.21. 
This corresponds to the slope index or 2/B in the notation of this paper. The mean break 
radius is 330 pc, which corresponds to only a few tens of millarcseconds at a typical lens 
redshift of z\ ^ 0.5. The bright images of strong lenses are therefore probing the steep outer 
part of the cuspy core profile. The light profile is much steeper than isothermal; however, 
the projected mass may be less steep depending on the distribution of dark matter. Central 
images are absent because break radii are small in cuspy core galaxies and so the steep slope 
continues to small radii. In fact, only 2 out of the 26 galaxies listed by Faber et al. (1997) 
have a sufficiently shallow profile for the central image to stand any chance of being visible at 
the 1% threshold. This provides an explanation of why most central images are unobservable 
at the current thresholds. 

Other than a central black hole, there is little evidence for dark matter inside the inner 
10 kpc of early- type galaxies that cannot be simply assigned to the stellar mass (e.g., Gerhard 
et al. 2001). Beyond 10 kpc or so, there probably is dark matter, although there are few 
hard facts on its distribution in ellipticals because of the absence of tracers at large radii. 
The dynamical evidence refers to the mass within spheres, whereas lensing is concerned with 
the mass within cylinders. Dark matter at large radii may alter the slope of the projected 
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mass distribution in the inner few kpc. However, provided the slope remains isothermal or 
steeper and the break radius remains small, then the central image remains unobservable. 

Fig. 11 shows the percentage of triplets with a visible central image as a function of 
the threshold for four different slope indices. Notice that for models with slopes steeper 
than isothermal, flattening typically makes the central image less visible, the reverse of the 
case when the slope is weak. The latter case shows the same trends as in the a — 0.01 
panel of Fig. 7. The radial caustic is large when B < 2 and the tangential caustic occupies 
a much smaller part of the space within it. Both fia and fii now decrease as q decreases, 
but Ha decreases faster to make the central image less visible. Fig. 11 shows us that the 
fraction of detectable images is a strong function of the threshold. It enables us to predict 
the threshold required to find the missing central images. Using the third panel of Fig. 11 as 
typical of the outer parts of cuspy core profiles, we see that the threshold has to be ~ 0.001% 
for the central image to be detectable in half the triplet systems. Even the most sensitive 
CLASS radio map (for B0218+357) probes only to a threshold of 0.05%, so this provides an 
explanation of why the central images have so far remained missing despite deeper searches. 



5. Conclusions 



We have shown how contour integration can ease the evaluation of the magnifications of 
images. This work extends the ideas presented in Hunter & Evans (2001) in two significant 
ways. First, our earlier analysis was restricted to scale- free power-law potentials. We have 
now generalized it to cover all elliptically stratified potentials. Second, we have obtained 
separate formulas for sums of the two direct images, for sums of the two inverted images, 
and for the magnification of the central image. Previously, we found only formulas for sums 
over all four images weighted with the signed magnifications. Our detailed applications are 
to power-law galaxies with cores. We have shown that the caustics are then always simple 
closed curves, and have given conditions for each of the four different caustic configurations. 
We have found an approximation for the magnitude of the central image which applies 
throughout the region inside the large outer radial caustic when the core is small and the 
slope index 2/B < 1. For small cores and weaker cusps, that approximation is directly 
useful only in a smaller inner region, but shows the magnitude of the central image to grow 
on a short length scale. We find that, for power-law lenses with small cores and an inner 
tangential caustic, the sums over separate pairs vary considerably with the image positions 
while the signed sums over all four images are generally remarkably uniform. That near- 
uniformity is a consequence of large cancellations between terms which vary with the position 
of the source. Hence uncertainties in the magnification can have major effects on estimates 
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of those four image sums. This lessens their usefulness for the modehng of lenses as proposed 
by Witt & Mao (2000) and as used also by us in our Paper I. We have shown that similar 
large cancellations can arise with any elliptically stratified potential, and not just power-laws. 
Hence, they will occur also with elliptically stratified densities in the limit of low eccentricity, 
and perhaps more generally, though this needs to be verified. 

As an application, wc have examined the constraints implied by the missing central 
images of triplet and quintuplet systems. There is only one convincing example of a grav- 
itational lens system with a central image, namely the triplet of the ultraluminous quasar 
APM08279+5255. However, there arc a total of ~ 50 doublets known. Although this is 
a heterogeneous sample, discovered by different observers using different techniques in dif- 
ferent wavebands, nonethless the probability of detecting central images does seem to be 
very low. We take a central image as observable if the magnification ratio of the faintest 
to the brightest image is ^ 1% (although some of the lens systems studied with the high- 
est sensitivity by CLASS do go deeper). A rough summary of the observations is that the 
probability of observing a central image for a triplet is ^ 2%. The absence of central images 
is understandable if the mass distribution in the lensing galaxy population is nearly cusped, 
and the cusp is isothermal or stronger. For typical source and lens redshifts, the size of the 
core radius a must be < 300 pc. The slope of the gravitational potential (3 must be < 1. 

Most of the optical depth to strong lensing resides in the most massive galaxies, namely 
giant ellipticals. Wc know from high resolution imaging of nearby giant ellipticals that they 
typically have cuspy cores with j3 ~ 0.7 or 2/B ^ 1.3 outside the break radius of a few 
hundred parsecs (Faber et al. 1997). The break radius corresponds to only a few tens of 
millarcseconds at a typical lens redshift. Hence, strong lensing is primarily probing the 
steep outer part of the cuspy core profile. This is much steeper than isothermal, as the 
surface density is falling typically like r~^'^. The cuspy cores by themselves can provide 
the explanation of the missing central images. Dark matter at large radii may alter the 
slope of the projected mass distribution in the inner few kpc. However, provided the slope 
remains isothermal or steeper and the break radius remains small, then the central image 
remains unobservable. The ratio of the faintest to the brightest image in cuspy core profiles is 
typically ~ 0.001%. Even the most sensitive radio maps available probe only to a threshold 
of 0.05%, so this explains why the central images have so far remained missing despite deeper 
searches. The sensitivity of the searches must be increased by a factor of ~ 50 to find them. 

NWE thanks the Royal Society for financial support. The work of CH is supported 
in part by NSF through grants DMS-9704615 and DMS-0104751. We thank an industrious 
referee for a helpful report. 
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A. Caustics 

In this Appendix, we show that the power-law galaxies with cores in the presence of 
external shear give rise to either one, three or five images, depending on the position of the 
source and the extent of the core radius. A point-like source is never lensed into more than 
five images. The caustics are simple closed curves and there is exactly one point on a caustic 
in each radial direction from the center of the lens. The models differ from elliptically 
stratified density distributions. For these, Witt & Mao (2000) and Keeton, Mao & Witt 
(2000) showed that external shear can cause butterfly and swallowtail cusps to develop on 
the caustics and so sextuple imaging and higher does occur. 

The caustics are curves in the source plane on which the imaging equation (23) has 
double roots. The partial derivative of the imaging equation with respect to t vanishes at a 
double root and hence we find caustics from the common solutions of the two equations 

t''[L{t)-a'K,{t)] = -K^{t), (Al) 
Bt''-\L{t) - a^Ki{t)] + t''[L\t) - a^K[{t)\ = -K[{t). (A2) 

Here, we have introduced the functions 

K^{t) = -{t - hf{t - hf, L{t) = [(Po - t)C - gC][(^o - t)C - QC\. (A3) 

Dividing the two sides of equations (Al) and (A2) gives a sextic polynomial in t, whose 
coefficients involve C,. Alternatively, we can solve for L{t) and L'{t) and divide the results to 
obtain a non-polynomial equation, but one which is independent of |C| and depends only on 
the angular argument (j) oi C,. This equation is 

K[{t) L'jty 

-2P^ + QP^\Q\ cos e-6P\Q\'^ + 2\Q\^ cos e l 

(P2-|g|2)(p2+|Q|2_2P|g|cos^) J ■ ^ ' 

As a — > 0, it reduces to equation (A2) of Paper I. It contains the complex angle 9 defined 
by QC/C — rather than 0. We find caustics by searching for roots for t along each 

angular direction in the source plane. Once t is known, |C| follows from equation (Al). In 
general, this search must be done numerically, but one can deduce from the graph of the 
term in square brackets (which has vertical asymptotes at t = ti and ^2 and is plotted in 
Figure 7 of Paper I) and the extra (1 — a^t^) factor, that there will always be a root in 
{t2,ti) for a point on the tangential caustic if a < 12^^"^ ■, and another in (ti,a~^/^) for a 
point on the radial caustic if a < There are four special directions for which points 



B 

1 



:i - aH""] 



- (l-a^O 
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on the caustics can be found explicitly. They are directions in which Qe is real, and we 
now consider them in turn. 

When Qe"^"^ = \Q\ and cos^ = 1, then t — t2 a>i which P = is a triple root 
of equation (A4). It gives a cusp at \C,\ — 2|(5|f^^'^^-\/l — a^ff on the tangential caustic 
provided a < ij^^^. The other solution of equation (A4) is the root of B{t—ti) — 2t{l — aH^) 
in (ti, a"^/"^) for a < t^^^^, which gives a point on the radial caustic. The t = t2 cusps on the 
tangential caustic are naked if their value of \(\ = 2\Q\t^^^'^ \J\ — a^t^ exceeds the value of 
1^1 for the fourth root for t of equation (A4). That value is given by equation (Al) evaluated 
for that root. The root cannot be found explicitly when a 7^ 0. Hence the marginal cases 
shown in Fig. 3 are found numerically by eliminating tjt\ between the equations 

B(t-ti) = 2t(l-aH^) 
t^{l-aH^){h-hf = t^il-aH^){t-hf, (A5) 

and solving for t2/^i for given B and a^tf . The second of equations (A5) comes from equating 
the two values for We can solve for t explicitly in the coreless a = case, to give the 
condition 



^ " ' ' , (A6) 



Bti 



for naked cusps. 



The other special directions are those for which Qe"^**^ = — and cos^ — —1. Then 
t — ti a,t which P — —\Q\ is a, triple root of (A4), and gives a cusp point on a caustic at 
\C\ — 2\Q\ti^^^^/l — aHf provided a < t'^^^'^. The other solution of equation (A4) satisfies 
B{t — 12) = 2t{l — aH^). Such a solution exists provided t2 < a~'^/^ . It gives a point on 
the radial or tangential caustic according to whether it is greater or less than ti. The t = ti 
cusp is on the tangential caustic in the first case, and on the radial caustic in the second, 
giving condition (21) as that for which the double lips configuration occurs. Double lips 
occur because the t = t\ cusp always lies on whichever is the inner caustic in this special 
direction; its |C| value never exceeds that for the other root of equation (A4). The two t 
and Id values are equal only when t — ti is a, quadruple root of equation (A4) and then 
B{ti — 12) — 2ti{l — aHf). This latter case is the transitional one in which the two caustics 
share a cusp and coincide locally, and marks the stage at which the region for core triplets 
has shrunk to zero. 

A simple analytical approximation for the radial caustic can be found for a small and 
B < 2. This caustic, which does not exist for a = 0, is large. It is found by looking for 
large roots of equations (Al) and (A2) for t, and approximating Ki by —t^ and L by i^|CP- 
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Eliminating t then gives the approximation 



B 



2-B 



{2-B)/B 



[l + 0(a2/^)]. 



(A7) 



Hence the radial caustic is approximately an ellipse elongated in the ?7-direction, as in Figs. 2b 
and d. The reason why those two ellipses for a = 0.1 are not well described by equation 
(A7) is that its limit of + = 1 is not accurate in the marginal case oi B = 2 until 
a -C 1. A more accurate formula for the radial caustic for i? = 2 is 



e + qV = 1 - 3a^/^ [(1 + 7i)f + 2^72^^ + q\l - liWf' + 0(a'/% 



(A8) 



We showed in Paper I, by an analysis of the quartic obtained when a = so there is no 
complicating power in equation (A4), that there is never more than one root in (^21^1) 
and one in {ti, 00). The additional (1 — a^t^) factor now present has only a small effect at 
finite t when a is small. However, it does ensure that equation (A4) always has a root in 
(ti,a~^/^) when this interval exists, which is not the case when a = if i? < 2. If more 
roots are to occur at larger values of a, there must be transitional cases at which the pair of 
equations (Al) and (A2) have a multiple root. We show next, using the sextic derived from 
those equations by eliminating the terms, that this cannot be, and hence that the cases 
listed in the previous paragraph are the only ones possible. 

To study the roots for t in t > ti, we work with the variable a — {t — ti)/\Q\, for which 
the sextic is: 

a^{a^ + 6(7^ + 12(7 + 8) + (B- 2)(7^ + 2[{B - 2)(2 + cos 6) 

-{po + cos9)]a^ + 2(1 + cose)[3{B-po - 3)a^ (A9) 
+2{B - 3po - 5)a - 4(1 +po)] = 0, 

where po = Po/\Q\ > 1. We let C„ denote the coefficient of cr" and apply Descartes' 
rule of signs (Henrici 1974) which shows that there is just one positive root for a when 
there is a single sign change in the coefficient sequence {Cq, C5, C4, C3, C2, Ci, Co)- We take 
(1 ± cos^^) > because the equations simplify in the special cases for which cos^^ = ±1 as 
discussed earlier and look for possible sign changes in the coefficient sequence. For B > 2, 
Cq, C5, and C4 are positive and Cq is negative. Also C2 > Ci as in Paper I, and C3 > if 
C2 > 0. For B — 2, Cq, C5, and C4 are positive, while C2, Ci, and Co are negative. For 
B <2, Cq and C5 are positive, C2, Ci, and Cq are negative, and C4 > C3. In each case, only 
one sign change can and does occur. 



a^BlQl 
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To show that there is just a single root in ^2 < ^ < ^i, we re-express the sextic in the 
variable t — {t — / {ti — t) which runs from to oo. It is then 

^ + (1 + 2r + r") (C4r^ + c^t^ + cir + cq) = 0, (AlO) 

where C4 > 0, C3 > 0, Ci, and Cq < are the coefficients of equation (A4) of Paper 1. Labelling 
this sextic as Cqt^ + C^t^ + C^t"^ + Cst^ + C2T^ + Cit + Co = 0, one finds that Ce, C5, and 
C4 are always positive and Co always negative. For B > 2, C2 > Ci and C3 > if C2 > 0. 
For B < 2, C2 and Ci are both negative. In either case, there is exactly one sign change in 
the coefficient sequence, and therefore exactly one root for t in (^2,^1)- 



B. Evaluation of the R Coefficients 

In this Appendix, we show how to evaluate the coefficients Ri{ti,t2,a^;j,i,k) and 
R2{ti,t2,a^;j,l,k) of the powers |A|^|i/|'^ introduced in Section 3. They can always be ex- 
pressed as finite sums. We consider first the scale-free case for which wc give compact explicit 
expressions, and then the cored case for which we derive recursive relations. We give explicit 
expressions for sums of coefficients Ri + R2 needed for four-image sums in the scale-free case. 
Then we derive a compact expression for the Ri + R2 sums valid for any elliptically stratified 
potential, and discuss its consequences. 



B.l. Scale- free Case 



Compact representations for the Ri and R2 terms can be derived for the coreless a — Q 
case by partial differentiation using Leibniz's rule. We find: 



i?i(ti,t2,0;j,£, k) 



1^ 



t 



Bj-1 



E 

m=0 



(h-t2)>'+^ 

jB - 1\ tf^-^-"^ 



\£— m— 1 



{e + k-m)\ 



m J{i-my. k\{ti - t2Y+''+^-"' 



This sum of {£ + 1) terms can be written as the terminating hypergeometric series 
i?i(ti,t2,0;j,£,A;) = l-^^-p^(^ ^ j,F,^-i,l-jB--k-i;l-^^ 



(Bl) 



(B2) 



with the proviso that, when A; = and the hypergeometric series becomes a geometric one, 
only the first {£ + 1) terms of that series are to be used. A similar expression for R2 can be 
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found from it by using the transformation relation in eq (36). It is 

^^2(^1, h, 0; J, £, k) = \ ^) 2F1 [-k, 1 - jB- -k-i-l-^j^. (B3) 

In obtaining separate formulae for the two pairs of bright images rather than for all four, 
the present work extends that of Paper I for the coreless a = case. We showed there that 
the coefficients for the sums of four images can all be expressed in terms of hypergeometric 
functions; that is we showed that 

-1 



Ri{h,t2,0]jJ,k) + R2{ti,t2,0;jJ,k) ^ 



(£ + k + l)\ 

k+e+l jB-k-l 



X 



n {jB-sf^^2Fi(^£+l,jB;£ + k + 2;l-f^. (B4) 



We also showed in Appendix C of Paper I that the infinite hypergeometric series (B4) can be 
represented as the sum of the two finite components (B2) and (B3). We did not appreciate 
the significance of the two separate components as two-image sums. We did warn of the 
tendency of the separate Ri and R2 components to cancel for small (1 — ^2/^1), and that 
their sum could be computed more easily using the rapidly convergent infinite series (B4). 
A simple instance of this is the i — case for which 



Ri{h,t2,0;j,0,k) 



{tl -t2)'+'' 



The sum here for R2 consists of the first {k + 1) terms in the infinite binomial expansion of 
{ti/t2)^^''^ in powers of (1 — ^1/^2); and hence the first {k + 1) terms in the expansion of 
—Ri. Another instance occurs when B — 1 or B — 2. Then Ri + R2 = for j > with the 
result that the sum of the four signed magnifications are then independent of the position of 
the source (Witt & Mao 2000; see also Paper I). These examples are simple instances of the 
large cancellations which can occur when contributions of the direct and the inverted image 
pairs, weighted with the signed magnifications, are combined. 



B.2. Cored Case 



The idea is to use Leibniz's rule to evaluate 



1 



Ri{ti,t2,a ■,3,i,k) = --^-Q^^ 



{h-t2f+^l-aHfy 



(B5) 
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The general derivative of the second factor is of the form 



fBj-l 



;i - aHfy 



,Bj-m-l 



(1 - aHfy 



n=0 



where the coefficients am,n, m > n > 0, can be found iteratively from the relation 

am,n = [B{j + n) - m]am-l,n + B{j + n- l)am-l,n-l, "0,0 = 1- 

Use of Leibniz's rule then gives 



Ri{ti,t2,a'^;j,i,k) 



e+i+Bj-i 



X 



E 



{i + k-my. ft2 



^ ml {i — my. \t 



n=0 



1-aHf 



In an analogous manner, we find; 
R2{ti,t2,a^;j,£,k) -- 



-lYt 



Bj-l 



- t2y+^+\i - aH^y 



X 



E 

m=0 



{i + k-my. fti 



m\{k — my. \t: 



- 1 



n=0 



The only a's for which it is simple to derive explicit expressions are 

B"(j + n-l)! 



an,o^Yl{Bj - s), 



OLr. 



(j - 1)! 



(B6) 



(B7) 



(B8) 



(B9) 



(BIO) 



The first set are the only ones that appear in the a = case when the double sums (B8) and 
(B9) reduce to the hypergeometric sums (B2) and (B3) respectively. 



B.3. CoefHcient sums 



We begin by rewriting equation (32) for R2 as 



R2{ti,t2,a'^;j,i,k) 



1 

mdt{dt\ 



G{t) (it 



2m Jc, {t-hyt-t2) 



(Bll) 



The choice G{t) — t^^ ^ / (1 — aH^y is needed for equation (32), but a similar equation with 
a different G will arise for some other eUiptically stratified potential ■0('?') with t — 2'0'(t). 



-42 - 



We again choose C2 to be a contour which encloses t — t2 but not t or any singularity of G. 
Then a residue calculation gives 

Git)dt Gjt^l^ 
ini L (t - hMt - 12) h-t2 ^ ' 



27ri Jc^ {t - ti){t - ^2) ^1 - ^2 ' 
There is a similar relation 

G{t) (it 



1 



2m Jc, {t-ti){t-t2) 



(B13) 



for the Ri coefficient. When we evaluate its integral by residues, and combine the results, 
we obtain 

^' + ^'- -mWM^^*'' ^ ■ ■ ■ ■ ■ ■ ^^^^^ 

Here 

G[t„t2] = ^^^i^^, (B15) 
ti — 12 

is a simple divided difference, while G[ti, . . . ti, ^2, ■ ■ ■ ^2] has ti repeated {I + 1) times and t2 
repeated {k-\-l) times, and is a divided difference of order + + 1) (Milne- Thomson 1960, 
chap. 1). 

The significance of formula (B14) is as follows. Divided differences are well-behaved 
when {ti — 12) becomes small. More specifically, a mean value theorem shows that 

R, + R^^ -G[h, ...h,t2,... t2\ = - ^^^j^^\y_ , (B16) 

where i, at which the {£ + k + l)th derivative of G is evaluated, is some value of t in the 
interval (^2,^1)- Thus, when the axis ratio q is close to 1 so that (ti — is small, while 
the function G varies smoothly in (^2,^1) as in Section B.l, then both Ri and R2 become 
very large when {ti — ^2) is small because both contain negative powers up to and including 
(^1 — ^2)"^"'^"^, while their sum varies little. But this is not the case with the lensing potential 
of Fig. 2a. Then a is large and the negative (1 — aH^y powers cause the derivatives of G to 
be large, so that i?i -|- R2 is then also large. 



